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We consider a reconstructing scheme using observational data from SNIa, BAO and CMB, based 
on a model of dark unification using a single non-minimally coupled scalar field. We investigate 
through a reconstruction program, the main features the current observational data imposes to the 
scalar field potential. We found that the form suggested by observations implies a step feature in 
the potential, where the kinetic and potential energy becomes of the same order of magnitude. 
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I. INTRODUCTION 

Dark energy is the name of the unknown component 
responsible for the current accelerated expansion of the 
universe In its simpler form, this can be described 
by a fluid with constant equation of state parameter w — 
— 1, corresponding to a cosmological constant, leading to 
the successful ACDM model, the simplest model that fits 
a varied set of observational data. 

This model posses a high dependence to initial condi- 
tions that makes it unnatural in many ways. For exam- 
ple, the current value for fl\ and VLm are of the same 
order of magnitude, a fact highly improbable, because 
the dark matter contribution decreases with a -3 , with 
a(t) the scale factor, meanwhile the cosmological con- 
stant contribution have had the same value always. This 
problem in particular is known as the cosmic coincidence 
problem. 

In this context, the most natural way to understand 
the acceleration of the universe, is to assume the exis- 
tence of a dynamical cosmological constant, or a theoret- 
ical model with a dynamical equation of state parameter 
(p/p = w(z)). The source of this dynamical dark energy 
could be both, a new field component filling the universe, 
as a quintessence scalar field [2J, or it can be produced 
by modifying gravity @. 

In [20j the authors suggested that the current obser- 
vational data favor a scenario in which the acceleration 
of the expansion has past a maximum value and is now 
decelerating. The key point in deriving this conclu- 
sion is the use of the Chevalier-Polarski-Linder (CPL) 
parametrization (2lj . p2| for a dynamical equation of 
state parameter 

w(a) = wo + (1 — a)wi, (I) 

where wo and w\ are constants to be fixed by observations 
and a is the scale factor. An update analysis performed 
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using recent SNIa data was informed in [25[ where similar 
conclusions were derived. Both analysis assume a flat 
universe. 

In [26| , we study the previous model, using a new set of 
data, the Union 2 data set [IH, and also considering the 
possibility to incorporate the curvature parameter fife , as 
a new free parameter in the analysis. We find that, the 
three observational test; SNIa, BAO and CMB, all can 
be accommodated in the same trend, assuming a very 
small value for the curvature, fife ~ —0.08. The best fit 
values suggest that the acceleration of the universe has 
already reached its maximum, and is currently moving 
towards a decelerating phase. 

Using a scalar field to model dark energy, it is possi- 
ble to reconstruct the scalar field potential from obser- 
vations. There are many approaches to do this Con- 
sidering models of a single non-minimally coupled scalar 
field, the rapid variations of the equation of state param- 
eter at low redshift can not be described. 

The quest for the unification of inflation, dark matter, 
and dark energy, in different combinations, by a single 
field, has been studied in Refs. The main mo- 

tivation behind all proposals is that we do not yet un- 
derstand the nature of the components responsible for 
the three phenomena, but we do know that their special 
properties are beyond the realm of the ordinary matter 
described by the Standard Model of Particle Physics. 

An extreme, most economical, possibility is that all 
three phenomena can be explained by the existence of 
one single field. As was first put forward in Ref. flH fioj) . 
the simplest option at hand is a scalar field <f> with a po- 
tential of the form V(4>) = V a + (l/2)m 2 2 . The energy 
scale Vq is to be set at the tiny value of the observed cos- 
mological constant considered in the concordance ACDM 
model, and the mass scale of the field, m ~ I0 -6 m p i, is 
determined by the amplitude of primordial perturbations 
generated during inflation. 

It is the purpose of this paper to explore the conse- 
quences of this trend, suggested from the observations, 
and look for special features through a characterization 
of a unification model based on a scalar field, that al- 
lows these transitions at low redshift. The reconstruc- 
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tion program uses as an intermediate phase a certain 
parametrization of w{z), which after its test using the 
data, is slightly altered looking for improvements in the 
fit. Here we use a \ 2 test considering the AIC and BIC 
criteria that controls the number of parameters to be 
used. 



II. w(z) FROM THE OBSERVATIONS 

The problem of extracting information of w(z) from 
observations can be understand in the following way. Be- 
cause most of the measurements give us information of 
the Hubble function H(z) or also the luminosity (or an- 
gular diameter) distance di,(z), we are forced to use 



w(z) = — — 



1 2(1 + z)HH' - 3H 2 



3i? 2 (l + z) 3 n m - H 2 ' 



(2) 



in the case of data from H(z), and we are forced to use 



w(z) 



2(1 + z)D" + 3D' 

d> - (D^n m (i + zy 



(3) 



in the case of the luminosity distance, where D = 
HqCLi,{z)/c. Notice that the precision in values for f2 m 
and Qk are crucial in this reconstruction procedure. 

The limitations of this process, first identified in 
are related to the dependence of the function w(z) of first 
and second derivatives of the H(z) and D(z) functions 
respectively. In order to tackle this problem, we can try 
to model the w(z) shape, by using a proper parametriza- 
tion. The most used is the already mentioned CPL [2lj |. 
[22l |. but there are many others designed to specific goals. 
For example, to describe a fast transition at rcdshift z t 
we can use 



w(z) =Wi + 



Wf 



l + cxp(^) 

or if we are interested in oscillatory behaviour we can use 

1 + Z C 

1 + z 



w(z) = wq + wi cos A lo. 



However, there is a concern regarding the number of pa- 
rameters used to parameterize w(z). It is clear that in- 
creasing the number of parameters [34j, is easiest to im- 
prove the fit with observations, but is not clear if in the 
meantime we are adjusting noise, instead of the truly 
physical relation. 



III. THE APPROACH 

In this section we describe the method to reconstruct 
the scalar field potential through the use of a iterative 
program improving the fit using two parameterizations 
for the equation of state parameter w{z). In this section 
we assume that fife = 0. 



The comoving distance from the observer to redshift z 
is given by 



-(*) = 



dz' 



where 



E 2 (z) 



Ho Jo E(z'Y 

n m {i + zf + n de f(z), 



(4) 
(5) 



exp < 3 



l + w(z') 
1 + z' 



dz' 



and fide = 1 — f^TO- The SNIa data give the luminosity 
distance di,(z) = (1 + z)r(z). We fit the SNIa with the 
cosmological model by minimizing the \ 2 value defined 
by 



XSNIa — ^ 



\p(zi) - Li bs{zi)} 2 



(6) 



fit 



where /i(z) = 5 log 10 [<iL(z)/Mpc] + 25 is the theoretical 
value of the distance modulus, and fi b s is the correspond- 
ing observed one. 

The BAO data considered in our analysis is the dis- 
tance ratio obtained at z — 0.20 and z — 0.35 from the 
joint analysis of the 2dF Galaxy Redsihft Survey and 
SDSS data 29], that can be expressed as 



£>v(0.35) 
£V(0.20) 



1.736 ±0.065, 



with 

Dv(zbao) = 



zbao 
H(zbao) 



dz 
H(zj 



2-1 1/3 



(7) 



(8) 



We fit the cosmological model minimizing the x 2 defined 
by 



2 _ [LV(0.35)/L>y(0.20) - 1.736] z 
Xbao — 



0.065 2 



(9) 



A result from the combination of SNIa and BAO is given 
by a joint analysis finding the best fit parameters that 



minimize Xsnu +X% AO- 



lux addition, we can incorporate to the analysis the 
CMB redshift parameter [3fJ, which is the reduce dis- 
tance at z ts = 1090 [H| 



R = Jn m H 2 r(zi s ) = 1.71 ± 0.019. 



We also apply the \ 2 



Xcmb 



[R- 1.71] 5 
0.019 2 



(10) 



(11) 



to find out the result from CMB and the constraints from 
SNIa+BAO+CMB are given by Xs N la + Xbao + Xcmb- 
Assuming a form for w(z) in terms of a certain number 
of parameters, we perform a Bayesian analysis to obtain 
the best fit values of all the free parameters in the model. 
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IV. THE CPL CASE 



For example, using the CPL parametrization, w = 
wq + wiz/(l + z), the function defined in (JS|) leads to 

/(*) = (1 + z f^o +W .) exp (-^j . (12) 

Using the Union 2 set |28j . consisting in 557 type la su- 
pernovae, the analysis leads to the results shown in Table 

in 



TABLE I: The best fit values for the free parameters using 
the Union 2 data set in the case of a flat universe model. 



Data Set 



X nun 



(I'd 



It' l 



SN 541.43 0.4197 -0.8632 -5.490 

SN+BAO 542.11 0.4281 -0.7959 -6.537 
SN+BAO+CMB 543.91 0.2547 -0.9979 0.190 



A plot of the deceleration parameter using these num- 
bers is shown in FigJTJ 



q(z) 




FIG. 1: Using the Union 2 data set we plot the deceleration 
parameter reconstructed using the best fit values for three 
cases: only SNIa (continuos line), SNIa+BAO (dashed line) 
and SNIa+BAO+CMB (dotted line) 




ion 



FIG. 2: The confidence limits for the case of Fig[T] 

The case SNIa alone and SNIa+BAO are almost iden- 
tical, but the case including CMB data does not show 



the rapid change at small redshift. Using the best fit val- 
ues for the parameters, give us the best fit function H(z) 
that we can use to reconstruct the scalar field potential. 
Using a standard procedure 33] with the equations for 
the flat case 



b\zf 
V(z) 



H' 



8irG 



(liz)M' 



(13) 
(14) 



we can plot directly the scalar field potential. In Fig[3]we 
show the integration of these relations. As is expected, 
in the cases SNIa and SNIa+BAO, we observe a rapid 
change in slope, a "knee" feature, that is not observed 
in the case including CMB. Then the sharp change ob- 




FIG. 3: Using the Union 2 data set we plot the scalar field po- 
tential reconstructed using the best fit function H(z) for the 
three cases: only SNIa (continuos line), SNIa+BAO (dashed 
line) and SNIa+BAO+CMB (dotted line). The first two over- 
lap almost exactly. 

served in the reconstructed deceleration parameter q(z) 
suggested by the data, is here visible in the scalar field 
potential V{4>) as this knee feature |35| . 

This is the anomalous effect mentioned first in [2(| us- 
ing the Constitution data set for SNIa, and also in [25[ 
and [26[ . A large negative value for w\ means that for 
small redshift the data suggest a negative slope for w(z). 
Because we are using the CPL parametrization, this fact 
spoils the large z behavior of the EoS parameter, leading 
to a large negative values w(z —> oo) ~ —7. 

This is the kind of small redshift transitions in w(z) 
that were discussed first in [23| and also in [24[. In (23j 
they use a form, 



w(z) = Wq + 



Wf — Wq 



l+CX P ((z-Zi)/A)' 



(15) 



that captures the essence of a single transition at z t in 
a range A, from wq initially to a final value Wf in the 
future. In [24| the authors discussed the possibility of a 
fast change in w{z) at z < 0.02 and its implication for 
a standard scalar field model. They found that while a 
canonical scalar field model can decrease the expansion 
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rate at low redshift, increasing the local expansion rate 
requires a non-canonical kinetic term for the scalar field. 

In this work we present an analysis of this problem 
using real data, as opposite to the previous analysis, to 
constraint the form of the equation of state parameter 
w(z) and then through relations (fT3|) to constraint the 
scalar field potential V(<j)) in a unified dark matter dark 
energy model. 



V. IMPROVING THE FIT 

As is evident from the previous section, the results 
show the incompatibility of the CPL parametrization in 
describing the variation of w(z) with redshift, because 
the data suggest a very large and negative value for w\ 
which spoils the large z behavior of w(z). 

In order to improve the fit, some authors have sug- 
gested the use of a new parametrization [20j] for w{z). 
We can use for example, the fast single transition form 
(|15| . to improve the fit. However, this parametrization 
has four parameters, two more than the CPL. Adding a 
new parameter would be justifiable only if the AIC or 
BIC (or a combination of these two) numbers indicate 
so [3J]. In this context, if we do not have new insight, 
we would like to keep the number of parameters fixed (in 
this case two), until one of these numbers (AIC or BIC) 
indicated something else. 

Let us start with a first iteration of the process. The 
negative value for w\ obtained in the previous section 
indicate that w(z) has to change its slope recently, and 
as the CPL reproduce perfectly, at large z its value does 
not change very much. So, let us use the following ansatz 
for w{z): 

Wo Wi 



w(z) = — 



l + z (1 + z) 2 



(16) 



This form has two parameters, like the CPL, so they 
are both directly comparable with % 2 statistic, and also 
this has a dust limit for z — > oo. The main feature of 
this parametrization is that allows the possibility to de- 
scribe a change of w(z) at low z, without spoiling the 
large z behavior. The result of a bayesian analysis is dis- 
played in Table|nl Notice that, although the deceleration 



TABLE II: The best fit values for the free parameters of the 
parametrization (|16|l using the Union 2 data set for SNIa, in 
a flat universe model. 



Data Set 


Xmin 




Wo 


Wl 


N 


541.32 


0.4121 


8.7568 


7.990 


SN+BAO 


542.11 


0.4281 


-0.7959 


-6.537 


SN+BAO+CMB 


543.91 


0.2547 


-0.9979 


0.190 



s 




FIG. 4: The same as FigfT] 




FIG. 5: The same as Fig[3] 



The second panel also shows this slight improve- 
ment, although still remains the conflict between low 
and high redshift data. The performance of this new 
parametrization is displayed in Table IIIII Clearly, our 



TABLE III: A summary of the performance using CPL and 
the new parametrization using only SNIa data. We display 
also the AIC and BIC indexes to compare with the ACDM 
fit. 



SNIa k -21nL AAIC ABIC 



ACDM 
CPL 

New 



542.63 
541.43 
541.32 




0.80 
0.69 




-1.2 
-1.31 



new parametrization is an improvement respect to CPL. 
The best fit parameters suggest a step like feature in the 
potential, as we can see in Fig. [6] We can still look for 
a better parametrization trying to accommodate all the 
three observational probes in a single trend. 



parameter obtained from SNIa and SNIa+BAO still re- 
mains very close each other, meanwhile that considering 
CMB behaves differently, the curve for the joint analysis 
SNIa+BAO+CMB shows a change in slope close to z = 
that is not observed using CPL. 



VI. RESULTS 

In this paper we investigate the consequences of the 
low redshift variations in the equation of state parame- 
ter, suggested by the data, and their implications in a 
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4> 

FIG. 6: The same as Fig. [3] Notice the improvement having 
all the curves overlap but still show a distinc behavior once 
we consider the CMB data. 

unified scalar field model for dark matter and dark en- 
ergy. This trend suggest a special feature in the scalar 



field potential. The reconstruction program uses as an 
intermediate phase a parametrization of w(z), which af- 
ter its test using the data, is slightly altered looking for 
improvements in the fit. Here we have used a \ 2 test 
considering the AIC and BIC criteria that controls the 
number of parameters to be used. We found a better 
parametrization than the CPL, that properly describe 
the low redshift behaviour, but although ameliorate the 
tension between low and high redshift, it can not describe 
all the data together. 
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